Electronic spectrum in high-temperature cuprate superconductors 
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A microscopic theory for electronic spectrum of the Cu02 plane within an effective p-d Hubbard 
model is proposed. Dyson equation for the single-electron Green function in terms of the Hubbard 
operators is derived which is solved self-consistently for the self-energy evaluated in the noncross- 
ing approximation. Electron scattering on spin fluctuations induced by kinematic interaction is 
described by a dynamical spin susceptibility with a continuous spectrum. Doping and temperature 
dependence of electron dispersions, spectral functions, the Fermi surface and the coupling constant 
A are studied in the hole doped case. At low doping, an arc-type Fermi surface and a pseudogap in 
the spectral function are observed. 
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I. INTRODUCTION 

Recent high-resolution angle-resolved photocmission 
spectroscopy (ARPES) studies revealed a complicated 
character of electronic structure and quasiparticlc (QP) 
spectra in copper oxide superconductors. In particular, 
a pseudogap in the electronic spectrum and an arc-type 
Fermi surface (FS) at low hole concentrations were re- 
vealed, a substantial wave-vector and energy dependent 
rcnormalization of the Fermi- velocity of QP ("kinks" in 
the dispersion) was observed (see, e.g. riband references 
therein). As was originally pointed out by Anderson 4 , 
strong electron correlations in cuprates play an essen- 
tial role in explaining their normal and superconducting 
properties. 

A conventional approach in describing strong electron 
correlations is based on consideration of the Hubbard 
models. The model has some advantages in comparison 
with the t-J model which can be derived from the Hub- 
bard model in the limit of strong correlations. Namely, 
the Hubbard model allows to study a moderate correla- 
tion limit observed experimentally in curates and more 
consistently takes into account a two-subband character 
of electronic structure, in particular, a weight transfer 
between subbands with doping. 

Various methods were proposed to study electronic 
structure within the Hubbard model. An unbiased 
method based on numerical simulations for finite clus- 
ters (for a review see e.g.|&) precludes, however, to study 
subtle features of QP spectra due to poor energy and 
wave-vector resolutions in small size clusters. In ana- 
lytical calculations of spectra mean-field type approxi- 
mations are often used (for a review seai£) which cannot 
reproduce the above mentioned effects caused by the self- 
energy contributions. In the dynamical mean field theory 
(DMFT) (for a review se o 9 i 10 ) the self-energy is treated 
in the single-site approximation which also unable to de- 
scribe wave-vector dependent phenomena. To overcome 
this flaw of DMFT, various types of the dynamical clus- 
ter theory were developed (for a review se o 11 ' 12 ). In these 
methods only a restricted wave-vector and energy reso- 



lutions can be achieved, depending on the size of the 
clusters, while the physical interpretation of the origin of 
an anomalous electronic structure in numerical methods 
is not straightforward. 

To elucidate the mechanism of the pseudogap for- 
mation, the charge carriers scattering by short-ranged 
(static) antifcrromagnctic (AF) spin fluctuations was 
considered in several analytical semi-phenomenological 
studies (for a review see£). More recently, by includ- 
ing into the DMFT scheme an additional momentum- 
dependent component of the self-energy originating 
from short-range AF (or charge) correlations, the spin- 
fluctuation scenario of the pseudogap formation 1 ^ and 
the arc-type FSi^ have been supported (for a review 
seei^). At the same time, it is important to study the 
effects of the charge carriers scattering by the dynami- 
cal spin-fluctuations which are believed to be responsible 
for the kink phenomenon^. This can be done by con- 
sidering the Dyson equation for the single-particle Green 
function (GF) within the Hubbard model in the limit 
of strong correlations. For instance, calculation of elec- 
tronic spectrum within the first order perturbation the- 
ory for the self-energy has reproduced quite accurately 
quantum Monte Carlo results^, while application of an 
incremental cluster expansion for the self-energy has en- 
abled to observe a kink structure in the QP spectrum— 

The aim of the present paper is to develop a micro- 
scopic theory for the electronic spectrum in strongly cor- 
related systems, as cuprates, which consistently takes 
into account effects of electron scattering by dynamical 
spin fluctuations. For this, we have considered an effec- 
tive Hubbard model reduced form the p-d model for the 
Cu02 plane in cuprates. By applying the Mori- type pro- 
jection technique for the thermodynamic GF— in terms 
of the Hubbard operators, we derived an exact Dyson 
equation, as was elaborated in our previous publica- 
tion o 19 i 20 i 21 . A self-consistent solution of the Dyson equa- 
tion with the self-energy evaluated in the noncrossing ap- 
proximation (NCA) beyond a perturbation approach was 
performed. 

This enabled us to calculate the dispersion and spec- 
tral functions of single-particle excitations, the FS, and 
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the electron occupation numbers. In particular, wc stud- 
ied a hole-doped case at various hole concentrations. At 
low doping, the FS reveals an arc-type shape with a pseu- 
dogap in the (7r,0) region of the Brillouin zone (BZ). A 
strong renormalization effects of the dispersion close to 
the Fermi energy ( "kinks" ) are observed due to electron 
scattering by dynamical AF spin fluctuations induced by 
kinematic interaction generic for the Hubbard operators. 
Electron occupation numbers show only a small drop at 
the Fermi energy. For high temperature or large hole con- 
centrations, AF correlations become weak and a crossover 
to a Fermi-liquid-like behavior is observed. 

In the next Section we briefly discuss the model and 
derivation of the Dyson equation, and the self-energy cal- 
culation in the NCA. The results of numerical solution 
of the self-consistent system of equations for various hole 
concentrations and discussion are presented in Sect. 3. 
Conclusion is given in Sect. 4. 



II. GENERAL FORMULATION 

A. Effective Hubbard model and Dyson equation 

Following a cell-cluster perturbation theory 
((\<^.^aa**i) based on a consideration of the orig- 
inal two-band p-d model for the Cu02 laye»24 we 
consider an effective two-dimensional Hubbard model 
for holes 
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where A" m = \in)(im\ are the Hubbard operators (HOs) 
for the four states n, m = |0), | ex) , 1 2) = | fj,), cr = 
±1/2 = (|, I), o = —a. Here s\ — Ed — fJ. and 
£2 = 2ei + U e ff where \x is the chemical potential. The 
effective Coulomb energy in the Hubbard model is 
the charge-transfer energy U e ff = A = e p — e^. The 
superscript 2 and I refers to the two-hole p-d singlet 
subband and the one-hole subband, respectively. Ac- 
cording to the cell-cluster perturbation theory, we can 
take similar values for the hopping parameters in (I): 
t" 
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„jj. The bare electron dispersion de- 
fined by the hopping parameter £y we determine by the 
conventional equation 



i(k) = 4i7(k) + 4i / 7 / (k) 



4t"j" 



(k), 



(2) 



where t, t', t" are the hopping parameters for the near- 
est neighbor (n.n.) (±a x , ±ay), the next nearest neighbor 
(n.n.n.) ±(a x ± a 9 ) and ±2a x ,±2a y sites, respectively, 
and 7(k) = (I/2)(cosfc 3J + cosfc y ), 7'(k) = cos k x cos k y 
and 7"(k) = (1/2) (cos 2k x + cos2fc y ) (the lattice con- 
stants a x = a y equal to unity). To get a physically rea- 
sonable value for the charge-transfer gap for the conven- 
tional value of t ~ 0.4 eV we take A = U e ff — 8t ~ 
3.2 cV. The bare bandwidth is W = 8t ~ U eff which 



shows that the effective p-d Hubbard model QJ corre- 
sponds to the strong correlation limit. In what follows, 
the energy will be measured in unit of t with e c i = 
in £1. The chemical potential pi depends on the average 
hole occupation number 
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The HOs entering obey the completeness relation 
Xf° + A,/ T + X^ + X? 2 = 1 which rigorously preserves 
the constraint of no double occupancy of any quantum 
state I in) at each lattice site i. Due to the projected 
character of the HOs, they have complicated commuta- 



tion relations 
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which results in the so-called kinematic interaction. The 
upper sign here refers for the Fermi-like HOs like Xf a 
and the lower sign is for the Bose-like ones, like the spin 
or number operators. 

To discuss the electronic structure within the 
model l|T]l. wc introduce a thermodynamic matrix Green 
function (GF)±& 



Gija (t — t ) 



((x ia (t)\xl(t>))) 

-i9(t-t')({X ia (t),xl(t')}}, (4) 



in terms of the two-component operators Xi a = 

Vcr2 \ 

X%a J and Xl = (Xf a Xf). To calculate the 

GF we apply the Mori-type projection technique by 
writing equations of motion for the Hcisenberg operators 
in the form: 



Zia — \Xi a , H] — EijaX^ 
j 



(5) 



where the irreducible Z-operator is determined by the 
orthogonality condition: 



{{zf;\ xD) = {z?;> x} a + x] a zr> = . (6) 



This defines the frequency matrix 

e ij = ({[X ia ,H],xl }) Q~\ 



(7) 



where Q = {{X la ,x} a }) = (^ Q Q 2 ° \ The weight fac- 
tors Q 2 = (Xf 2 +XD = n/2 and Q 1 = (X™ + Xf°) = 
1 — Q2 in a paramagnetic state depend only on the hole 
occupation number The frequency matrix (JJJ deter- 
mines the QP spectra within the generalized mean field 
approximation (MFA). The corresponding zero-order GF 
in MFA reads: 



d°(k,u) = (wfo-e(k)) l Q, 



(8) 



where fo is the unity matrix and we introduced the fre- 
quency matrix (JJJ in the k- representation e(k). By dif- 
ferentiating the many-particle GF ((Zl a r (t) | X f ja (t'))) 
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over the second time t' and applying the same projection 
procedure as in JSJ we derive the Dyson equation in the 
formic 



G^k,^)" 1 = G^k.w)" 1 - E ff (k,w) 



(9) 



Here the self-energy matrix E CT (k, ui) is determined by a 
proper part (which have no single zero-order GF) of the 
many-particle GF in the form 

s CT (k,u/) = crWW^LT ( 10 ) 

The equations JHJ ~ l|10|l provide an exact representation 
for the GF J3J. However, to calculate it one has to use an 
approximation for the self-energy matrix 1)10(1 which de- 
scribes inelastic scattering of electrons on spin and charge 
fluctuations. 

It is important to point out that in the Hubbard 
model Q electron interaction with spin- or charge fluc- 
tuations are induced by the kinematic interaction with 
the coupling constants equal to the original hopping pa- 
rameters, as has been already pointed out by Hubbard 5 . 
For instance, the equation of motion for the operator Xf 2 
reads 

idXf/dt = \Xf,H\ = {e x +A)X? 2 



E 



1 ,22 R 22 
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where B°z , are Bose-like operators describing the num- 
ber (charge) and spin fluctuations: 

Biaa> = (Xf 2 + X^ cr )5 l ji a + X° a 5 a , s 

= (N l /2 + S~)5 rr , a + S°6 a ^, (12) 

Therefore, in the Hubbard model QJ, contrary to spin- 
fermion models where electron interaction with spin- or 
charge fluctuations are specified by fitting coupling con- 
stants*, this interaction is fixed by the hopping parame- 
ters. 

B. Mean-Field Approximation 

The single-particle excitations in MFA are defined by 
the frequency matrix J7J). By using equations of motion 
like {HJ, we get the following energy spectrum for holes 
in two subbands 

ei,a(k) = (l/2)Mk)+ Wl (k)] T (l/2)A(k), 
A(k) = {[ W2 (k)- Wl (k)] 2 +4^(k) 2 } 1 /2 j (13) 

where the original excitation spectra in the Hubbard sub- 
bands and the hybridization parameter are 

wi(k) = 4iai7(k) + 4i / /3i7 / (k)-M, 
w 2 (k) = 4ta 2 7(k)+4t'/3 27 '(k) + A-^, 
W(k) = 4tai27(k)+4i'/? 12 7'(k). (14) 



where we omitted t" contribution in J2J) and in- 
troduced the renormalization parameters £*i(2) = 

Ql(2)[l + Cl/Qi(2)]' A(2) = Ql(2)[l + C2/Q 2 (2 )] , "12 = 

y/QZUl-Ch/QM 0i2 = VQIQ^[1-C 2 /Q 1 Q 2 }. As 
in the Hubbard I approximation, we neglect number fluc- 
tuations (dNiSNj)^j-j but take into account contribu- 
tions from the spin correlation functions for the n.n. and 
the n.n.n. sites: 



Gi = (SiS 



i±a x /a B / ) 



Co 



(SiS 



(15) 



The renormalization of the QP spectra l|13l) , (|14|) caused 
by strong spin correlations in the underdoped region re- 
sults in suppression of the n.n. hopping which changes 
the shape of the spectra and reduces the bandwidth. For 
instance, if we consider the limiting case of the long- 
range AF Neel state with the n.n. correlation function 
Gi ~ -1/4 at half-filling, Qi = Q 2 = 1/2, we obtain 
a i(2) = 0. This results in complete suppression of the 
n.n. hopping and transformation of the spectra HI 411 to 
the n.n.n. hopping oc t'7'(k) as was discussed iniS. 

For the diagonal components of the zero-order GF (JSJ 
we have 



c ; la|(u ,, J-»''-;^ WW, (16) 

( ' w-ei (a) (k) w-£ 2(1) (k)' 



where the parameter 

, nA _ gg(k) - cjg(k) _ 1 cj 2 (k) -o;i(k) 

6(k, ~ £2 (k)- £l (k) ~2 ~~ 2ACk) ( 7) 

determines the contribution due to the hybridization. 

C. Self-energy Corrections 

Dyson equation (O for the GF is convenient to write 
in the form 

G ff (k,u) = (a;fo-£(k)-S„(k ) a))) * Q, (18) 

where the self-energy reads 

S CT (k, W ) = ((if)|^f)t))^)Q- 1 . (19) 

To make the problem tractable, we can neglect in 
the self-energy matrix l|19|) the off-diagonal components 
Si2,o-(k, to) in comparison with the hybridization parame- 
ters W(k) in l|14|). This enables us to write the diagonal 
components of the full GF (|18fl in the form similar to 
GEJ: 



Gn(22)(k, LU) = 



Qi( 2 ) [1 - 6(k)] 



ui - £i( 2 )(k) - S U ( 22 )(k,a;) 

Qi (2 ) b(k) 
> ~ £2(1) (k) - £22(11) (k w) 



(20) 
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Here the hybridization parameters b(k) are determined 
by the formula similar to (|17f) which gives an accurate 
approximation for a small doping at n ~ 1 . 

Now we calculate the self-energy (|19|l in the non- 
crossing (NCA) or the self-consistent Born approxima- 
tion (SCBA) by neglecting vertex renormalization. As 
follows from the equation of motion the op- 

erators determined by J5J) are essentially a product of a 
Fermi-like Xj(t) and Bose-like Bi(t) operators. In SCBA, 
the propagation of these excitations of different types in 
the many-particle GF in (|19|l are assumed to be indepen- 
dent of each other. Therefore, they can be decoupled in 
the time-dependent correlation functions for lattice sites 
{i 7^ h I 7^ m ) as follows 

(fli(t)X,-(t)B,X m ) ~ {X j {t)X m ){B i {t)B l ). (21) 

Using the spectral representation for these correlation 
functions, we obtain the following formula for the di- 
agonal self-energy components Su(22)(k,w) = S(k, ui) 
which are the same for two subbands: 

£(k,w) = jjYl /^(w,«|q,k-q) 
q 

x (-l/7r)Im[Gi(q,z)+G 2 (q,z)], (22) 
where the corresponding subband GFs are: 

1 



Gi ( 2)(q,w) 



w - £i( 2 )(q) - S(q,w) 



(23) 



The kernel of the integral equation i12L'li has the following 
form: 



K(a J ,z\ q ,k-q) = \t( q )f J 



an 



uj — z — Q 



x [tanh(^/2T) + coth(fi/2T)]Imx sc (k-q,fi), (24) 

where the interaction is defined by the hopping parame- 
ter t(q) (J2J). The spectral density of bosonic excitations is 
determined by the dynamic susceptibility of the Bose-like 
operators Bi(t) in (|21() - the spin and number (charge) 
fluctuations: 

XsM^) = -[«S q |S_ q » w + (1/4){(6N«\6N-M , (25) 

where we introduced the commutator GF for the spin S q 
and the number SNq = — (AT q ) operators. 

Thus we obtain a self-consistent system of equations for 
the GFs l|23[) an d the self-energy l|22|) . A similar system 
of equations was obtained within the composite operator 
method^. In comparison with the t-J model studied by 
us inSH, for the Hubbard model JU we have two contribu- 
tions in the self-energy l122[| determined by the two Hub- 
bard subbands, while in the t-J model only one subband 
is considered. However, depending on the position of the 
chemical potential, a substantial contribution to the self- 
energy comes only from the GF of those subband which 



is close to the Fermi energy. A contribution from the GF 
of another subband which is far from the Fermi energy, is 
suppressed due to a large charge-transfer energy A in the 
denominator of those GF. Neglecting the latter contribu- 
tion, we obtain a self-consistent system of equations for 
one GF close to the Fermi energy and the corresponding 
self-energy function similar to the t-J model—. 



III. RESULTS AND DISCUSSION 

A. Self-consistent system of equations 

To solve the system of equations for the self-energy 
(12211 and the GFs 12-ill we should specify a model for the 
spin-charge susceptibility (|25() . Below we take into ac- 
count only the spin- fluctuation contribution Xs(q, = 
— {{S q | S_ g )) w for which we adopt a model suggested in 
numerical studies^ 



Imx.s(q,w + i0 + ) = x«(q) X s H 



Xn 



i + e 2 (i + 7(q)) 



tanh 



1 



2Tl + (uj/uj s y 



(26) 



The q-dependence in Xs (q) is determined by the AF cor- 
relation length £ which doping dependence is defined be- 
low. The static susceptibility xo at the AF wave vector 
Q = (7T, 7r) is fixed by the normalization condition: 



1 

7T 



+oc 



— /An r xl'(z) 4y> 8 (q), (27) 

exp (z/T) — 1 N 



which gives the following value for this constant: 



Xo 



o; s (S?) {iV^l + ^[i + 7 ( q )]} 



(28) 



In H3) we introduced (Sf) = 3(5| Sf) = (3/4)((l-X°°- 
X? 2 )) ~ (3/4)(l - \S\) where at the hole doping 5 ~ 
{X? 2 ), while at the electron doping 6 ~ —(Xf°). 

The spin correlation functions i|15|l in the single- 
particle excitation spectra (|13|l in MFA are defined by 
equations 

C i = j?}2 C «^ ^ = ^£C q7 '(q). (29) 
q q 

The static correlation function C q can be calculated 
from the same model ()26[) as follows 

c(0 



C a 



(SqS- q ) 



i + e 2 [i + 7(q)] ' 

where the factor C(£) = xo (w a /2). 



(30) 
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TABLE I: Static spin correlation functions JSS), C(£) QOt 
and the AF correlation length £ in Il2(it at various hole con- 
centrations n = 1 + 5 



5 = 



0.03 0.05 0.10 0.15 0.20 0.30 



Ci -0.36 -0.26 -0.21 -0.18 -0.14 -0.10 

C" 2 0.27 0.16 0.11 0.09 0.06 0.04 

(7(C) 22.0 5.91 3.58 2.67 1.93 1.40 

£ 8.0 3.40 2.50 2.10 1.70 1.40 



To specify the doping dependence of the AF correlation 
length £(<5) at low temperature, we fit the correlation 
function C\ calculated from (|29|l to the numerical results 
of an exact diagonalization for finite clusters — . The 
values of the AF correlation length, calculated values of 
C 2 and the correlation function C(£) = (S q S_ q ) at the 
AF wave- vector q = Q = (it, ir) are given in Table |I] 

To perform numerical calculations, we introduce the 
imaginary frequency representation for the GF (j23|) : 



Gl(2)(q, iu> n ) 



£i(2)(q) - S(q, iu)r, 



(31) 



where iuj n = inT(2n+ 1), n = 0, ±1,±2, ... . For the 
self-energy (|22[1 we obtain the following representation: 



S(k,iw n ) = -— y^y^[6?i(q,io; m ) +(7 2 (q,ta; m )] 

q m 

x A(q, k - q | iu> n - iiv m ) . (32) 



The interaction function is given here by the equation 
A(q, k - q | iu v ) = -\t(q)\ 2 Xs (k - q) F s (iu> u ), (33) 



where the spectral function: 



1 



2xdx 



x 2 + {uo v /uj s ) 2 1 + X 2 



tanh 



XLO s 

~2T 



(34) 

Let us compare the self-consistent system of equations for 
the GF (|3"TJl and the self-energy l(3"2"|) with results of other 
theoretical approaches. In our theory based on the HO 
technique we start from the two-subband representation 
for the GF Q which rigorously takes into account strong 
electron correlations determined by the Coulomb energy 
U e ff. This results in the Mott gap at large U e tf (sec 
below) as in the DMFT. On the other hand, the kine- 
matic interaction, generic to HOs, induces the electron 
scattering by spin (charge) dynamical fluctuations l|25(l 
which are responsible for the pseudogap formation as in 



12.27 



or 



the two-particle self-consistent approach (TPSC) 
the model of short-range static spin (charge) fluctuations 
- the £k-modei£. 

To prove this, let us consider the classical limit for the 
self-energy (I32|l by taking into account only the zero Mat- 
subara frequency iuj v = in the interaction (|33|) which 



gives ioj m = ioj n in i|32|) . In the limit of large AF corre- 
lation length £ 3> 1 the static spin susceptibility x s (q) 
in H26fl shows a sharp peak close to the AF wave-vector 
Q = (jr, 7r) and can be expanded over the small wave- 
vector p = q Q: 



Xs(q) 



\0 



ep 2 



(35) 



where we introduced K = £ 1 and took into account that 
the constant igSJ xo - A£ 2 with A = (87r/w s )(S?)[ln(l+ 
47r £ 2 )] -1 for the square lattice. In this limit we get the 
following equation for the self-energy l(3"2"|) : 

£(k^ n )^| g (k-Q)| 2 ^]T^^ 

p 

x [Gi(k - Q - p, iivn) + G 2 (k Q p, iw„)], (36) 



where the effective interaction 



(37) 



|. 9 (q)| 2 = yl|t(q)| 2 ^(0). 

Expanding the QP energy £ 1 ( 2 )(k— Q — p) ~ 
£i(2)(k — Q) — p-v 1 ( 2 \ k _Q we obtain for the GFs 
in 13ti|) the following representation: 

Gi( 2 ) (k Q p, ioj n ) ~ {iu n ~ £i( 2 ) (k - Q) 
+ p-v 1(2) , k _ Q -S(k-Q,^„)} _1 - (38) 

The system of equations for the GFs (|38J) and the 
self-energy (|36|l is similar to those one derived in the 
TPSC approach 2 ^) and the Sk-model^ apart from the 
interaction function and the two-subband system of 
equations. In our approach the vertex l|37|l is deter- 
mined by the hopping parameter |t(k — Q)| 2 , while in 
the TPSC and the Sk-model the coupling constant is 
induced by the Coulomb scattering, e.g., ini^ g 2 = 
U 2 ((niiriii)/n 2 )(S 2 ) /3. However, the values of these 
ver tices a re close: the averaged over the BZ value 
(-\/|i(k)| 2 )k ~ 2t is comparable with the coupling con- 
stant g < 2t used ini^. In the spin-fermion model the 
self-energy is also determined by spin-fluctuations (see, 
e.g.^) with the coupling constant fitted from ARPES 
experiments g ~ 0.7 cV^ 2t of the same order. As 
in the TPSC theory, in the limit £ — > oo the AF gap 
A J 4j?(k) oc |t(k — Q)| 2 in the QP spectra emerges in the 
subband located at the Fermi energy. This result readily 
follows from the self-consistent equations for the GF l|31(l 
with the self-energy (|36|l where in the right-had side GF 
(|38|l is taken at p = 0. Thus, in our approach the pseudo- 
gap formation is mediated by the AF short-range order 
similar to TPSC theory and the model of short-range 
static spin fluctuations in the generalized DMFT^. 

In the next sections we consider the results of self- 
consistent calculations of the GFs H31|) and the self- 
energy (|32|l in the hole doped case for various hole con- 
centration 5 = n - 1 > 0. In Sects. IhTbI - MI PI the 
calculations were performed at temperature T = 0.03t ~ 
140 K and T = 0.3t for A = 8t, t ~ 0.4 cV and 
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t' = — 0.3t. Several results are reported for A = At, t' = 
-0.13i, t" = 0.16* in Sect. EH! For the spin-fluctuation 
energy in (j26|) we take co s = OAt. The AF correlation 
length and the static correlation functions C\, C2 in 
(|T5|) are defined in Table U 



B. Dispersion and spectral functions 

In ARPES measurements and QMC simulations the 
spectrum of single-electron excitations is determined by 
the spectral function ^4( e z) (k, uS) = -A(h) (k, —oS) . The 
spectral function for holes can be written as follow: 

Aft)( k > w ) = --Im((aka |oL))^+iO+ 
= [Qi + P(k)]Ai(k,w) + [Q a - P(kp 2 (k,u;).(39) 

Here we introduced for the hole annihilation ak CT and cre- 
ation operators the definition in terms of the Hub- 
bard operators a ka = Xf a + 2a Xf 2 , aj^ = Xf° + 
2aX 2a and used all four components of the matrix GF 
(|T£|l G a (i (k, lS) with the diagonal components given by 
(j20(l . In l|39|l we introduced also the one-band spectral 
functions determined by the GFs l|23|) : A 1 ( 2 )(k, cj) = 
— (I/71") ImG 1 ( 2 )(q,w) • The hybridization effects are al- 
lowed for by the parameter P(k) = (?i — l)6(k) — 
2VQl^W(k)/A(k). 

The dispersion curves given by maxima of spectral 
functions l|39|) were calculated for hole doping S = 0.05 — 
0.3. At low hole doping, 5 = 0.05, 0.1, the dispersion 
reveal a rather flat hole-doped band at the Fermi energy 
(FE) (oj = 0) as shown in the upper panel in Fig. ^ 
The corresponding spectral function (the bottom panel) 
demonstrates weak QP peaks at the Fermi energy. With 
doping, the dispersion and the intensity of the QP peaks 
at the Fermi energy substantially increase as demon- 
strated in Fig.Hthough a flat band in X(ir, 0) -» T(0, 0) 
direction is still observed in accordance with ARPES 
measurements in the overdoped Lai.7sSro.22CuO^£. To 
study an influence of AF spin-correlations on the spectra, 
we calculate the spectral functions at high temperature 
T = Q.3t for 8 = 0.1 by neglecting spin correlation func- 
tions (|15|) in the single-particle excitation spectra 113fl in 
MFA and taking a small AF correlation length (£ = 1.0) 
in the spin-susceptibility 1|26|) . Figure[3]shows a strong in- 
crease of the dispersion and the intensity of the QP peaks 
at the Fermi energy as in the overdoped region, 8 = 0.3, 
which proves a strong influence of AF spin-correlations 
on the spectra. A crude estimation of the Fermi veloc- 
ity from the dispersion curve in the r(0,0) — * M(tt, it) 
direction in Fig. [3 for the overdoped case gives the value 
Vf — 7.5t A ~ 3 (eV-A) for the hopping parameter 
t = 0.4 eV which can be compared with experimental re- 
sults V F 21 2.2 (eV-A) for overdoped Lai. 78 Sr .22CuO/r& 
and V F ^ 3.9 (eV-A) for overdoped Bi-221222. With 
doping, the electronic density of states (DOS) shows a 
weight transfer from the upper one-hole subband to the 
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FIG. 1: Dispersion curves (upper panel) and spectral func- 
tions (bottom panel) in units of t along the symmetry direc- 
tions T(0, 0) -> M(tt, tt) -» X(ir, 0) -> T(0, 0) for 8 = 0.05. 
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FIG. 2: The same as Fig. Q for hole concentration 8 — 0.3. 



lower two-hole singlet subband as shown in Fig. How- 
ever, even in the overdoped case a noticeable part of the 
DOS retains in the upper one-hole subband. 

It is interesting to compare our results with those ob- 
tained in the generalized DMFT— which should be close 
to each other as discussed at the end of Sect. IIII Al 
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cluster perturbation theory^ (see Fig. 2 (a) ii 
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C. Fermi surface and occupation numbers 

The Fermi surface for the two-hole subband was deter- 
mined by a conventional equation: 



e 2 (k F ) + Re£(k F ,w = 0) = 0, 



(40) 



as shown in Fig. [5] and then compared with those 
one obtained from maxima of the spectral function 
A e i(k 7 uo = 0) on the (k x , fc. y )-plane for S = 0.1,0.2 
shown in Fig. [S] The FS changes from a hole arc-type 
at 6 = 0.1 to an electron-like one at 8 = 0.3. Exper- 
imentally an electron-like FS was observed in the over- 
doped Lai.78Sro.22CuO,f2&. The doping dependent FS 
transformation can be also observed by studying the 
electron occupation numbers. The electron occupa- 



FIG. 3: The same as Fig. but for the hole concentration 
S = 0.1 and at high temperature T = 0.3t. 




Energy 

FIG. 4: (Color online) Doping dependence of the electronic 
density of states. 



In fact, the spectral function shown in Fig. 8 for 
t' = —0.4 demonstrates a similar flat QP bands in 
r(0, 0) -> X(tt, 0) and r(0, 0) ->■ M (tt, tt) directions, as in 
our Fig. ^ and Fig. [21 a strong intensity transfer from the 
lower electronic Hubbard band (LHB) to the upper Hub- 
bard band (UHB) at the M(w, tt) point of the BZ and a 
splitting of the LHB close to the X(ir, 0) point. An analo- 
gous temperature and doping (£) behavior of the spectral 
functions and the pscudogap revealed in the both theories 
supports the spin-fluctuation scenario of the pseudogap 
formation. A similar behavior was observed also in the 
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FIG. 5: (Color online) Doping dependence of the FS for S = 
0.1 (full line at T = 0.0M and dotted line at T = 0.3i), 8 = 0.2 
(dashed line), and 5 = 0.3 (dot-dashed line). 





FIG. 6: A(k,a; = 0) on the FS for S = 0.1 (left panel), 5 = 0.2 
(right panel). 

tion numbers in (k)-space for one spin-direction equal 
to Af( e ;) (ct, k) = 1 — Nff l ){ (J i k) where the hole occupation 
numbers JV^^ct, k) = Nt^Qt) according to Q are deter- 
mined only by the diagonal GFs (|20|l . From the latter 
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equation and i|23|) we get: 



N {h) (k) = [Q! + (n-l)6(k)]JVi(k) 
+ [Q 2 -(n-l)6(k)JV 2 (k), 




FIG. 7: (Color online) The electronic occupation numbers 
iV k for <5 = 0.1 at T = 0.03t (upper panel) and at T = 0M 
(bottom panel). 




FIG. 8: (Color online) The electronic occupation numbers A?k 
at T = 0.03i for 5 = 0.3. 



The electron occupation numbers in a quarter of the 
BZ (0 < k x ,k y < 7r) are shown in Fig.[7]for 5 = 0.1 at 
low temperature T = 0.03i and at high temperature T = 
0.3t. With doping the the shape of the iVk is changing 
revealing a transition of the holedike FS to the electron- 
like in the overdoped case 8 = 0.3 as plotted in Fig. EI 

While in the underdoped case at 6 = 0.1 the drop of 
the occupation numbers at the Fermi level crossing is 
rather small, ANun — 0.15, for high temperature T = 
0.3i or in the overdoped case at S = 0.3 when the AF 
spin correlations are suppressed, the occupation number 
drops are substantially increased: AN^ ~ 0.45, 0.55, 
respectively. Thus, the arc formation and a small change 
of the electron occupation numbers at the FS crossing at 
low doping further prove a large contribution of the spin 
correlations in the renormalization of QP spectra. 

The obtained result concerning the "destruction" of 
the FS caused by the arc formation shown in Fig. [BJ and 
Fig. 1161 for low doping, which corresponds to large £, 
correlates well with the studies within the generalized 
DMFT— . As shown in Fig. 2 inii, the spectral density 
intensity plots clearly demonstrate the arc formation on 
the FS for large coupling constant A s / = A = 2t and 
£ = 10, while the FS determined from (|40|l gives several 
solutions as in our Fig. El for U e ff = it in Scct. lIIIEl 



D. Self-energy and kinks 

Energy dependence of the real and imaginary parts of 
the self-energy S(k,w) for 6 = 0.1, 0.3 at the r(0,0), 
S(tt/2,it/2) and M(tt, 7t) points are shown in Fig. 
These plots demonstrate a strong dependence of the 
self-energy on the wave-vector and the hole concentra- 
tions. With doping, the coupling constant substan- 
tially decreases as seen by the decreasing of the imag- 
inary part and the slope of the real part at the FS 
crossing which determines the coupling constant A = 
— (<9Re£(k, Lj)/<9tj) w= o . As shown in Fig. 1101 the cou- 
pling constant in the r(0, 0) — > M(ir, it) direction de- 
creases from A ~ 7.86 at 5 = 0.1 to A ~ 3.3 at S = 0.3. 

At large binding energies (greater than the boson en- 
ergy responsible for the interaction) the self-energy ef- 
fects vanish and the electron dispersion should return to 
the bare value, giving a sharp bend, the so-called "kink" 
in the electron dispersion. The amplitude of the kink 
and the energy scale where it occurs are related to the 
strength of the electron-boson interaction and the boson 
energy, respectively. In ARPES the kink is observed as a 
changing of the slope for an intensity plot for the spectral 
function A(k, u>) in a particular k-wave vector direction 
below the Fermi level u> < (for electrons). Usually 
two directions are studied: the nodal (r — ► M) and the 
antinodal (X — > M) ones. Intensity plots for the spectral 
function A(k,uj) at S = 0.1 are shown in Fig. ^] in the 
nodal direction (left panel) and the antinodal one (right 
panel). The same plots at 8 = 0.3 are shown in Fig. 1121 
in the nodal direction (left panel) and X(ir, 0) — > T(0, 0) 
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FIG. 9: (Color online) Energy dependence of the real and 
imaginary parts of the self-energy E(k, w) at the F(0,0), 
S(ir/2, 7r/2) and A/(7r,7r) points at 8 — 0.1 (upper panel) 
and S = 0.3 (bottom panel). 



direction (right panel) . A change of dispersion is clearly 
seen with increasing binding energy below the FS shown 
by dotted line. For the underdoped case the kink is larger 
than for the overdoped one. A crude estimation of the 
strength of the kink from the ratio of the dispersion slope 
Vf close to the FS (u) = 0) to those one Vp at large bind- 
ing energy (u) ~ 0.2t), Vg/V F = (1 + A), gives the fol- 
lowing values: (1 + A) ~ 7.6, 3.5 at S = 0.1 for the nodal 
and antinodal directions, respectively. In the overdoped 
case the nodal value is much smaller, while in the the 
antinodal X(ir,0) — > r(0,0) direction is still quite large: 
(1 + A) ~ 2.5. These estimations are in accord with the 
evaluation of the coupling constant A from the slope of 
the real part of the self-energy discussed above. 

It is important to stress that in our theory the self- 
energy effects and the corresponding kinks are induced 
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FIG. 10: (Color online) Re£(k,w) in the r(0,0) -> M{tt,it) 
direction at the FS. 
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FIG. 11: (Color online) Dispersion curves along the symme- 
try directions M(tt, it) — * r(0, 0) (left panel) and M(ir, n) — > 
X(ir,0) (right panel) in units of t for 5 = 0.1, T = 0.03t. 
Fermi level crossing is shown by vertical dotted line. 
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FIG. 12: (Color online) The same as Fig.[TT]but for S = 0.3 
along the symmetry directions M(n,n) — > r(0,0)(left panel) 
and X(iv,0) T(0, 0) (right panel). 



by the spin-fluctuation spectrum in the form of the con- 
tinuum l|26|) which at low temperature T ~ 0.03t <C 
uj s = OAt has a large intensity already at small en- 
ergy to ~ 0.03t and decreases slowly up to a high energy 
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w ~ t. In the spin-fcrmion model the kink phenomenon 
is usually explained by electron interaction with the spin- 
resonance mode f2 ros — 40 meV observed in the supercon- 
ducting state. This results in a break of the electron dis- 
persion (" kink" ) at a certain energy oj ~ f2 res + Ao where 
Ao is the superconducting gap (see, e.g£). In the normal 
state considered in our theory the spin-resonance mode 
is inessential. Its contribution amounting only few per- 
cents of the total spin fluctuation spectrum (|27() should 
not change our results which reveal a rather strong in- 
teraction with a smooth energy variation without any 
specific kink energy. 



E. Dispersion and FS at U eH - = A = 4t 

The effective Coulomb energy in the Hubbard 
model JU U e ff = 8t results in a large charge-transfer 
gap A ~ 3 eV for t = 0.4 eV even in the overdoped case, 
Fig. |2 while experiments point to a smaller value of the 
order of 1.5 — 2 eV. To correct this inconsistency, we 
present in this section the results obtained for a smaller 
value of U e ff = A = At. We also take into account the 
hoping parameter for the n.n.n. ±2a x , ±2a y sites and 
fix the hoping parameter in the model dispersion J2J) as 
suggested for the effective Hubbard model based on the 
tight-binding fitting the LDA calculations for L^CuO^Si 
as follow: t' = -0.13*, *" = 0.16* with * ~ 0.7 eV. 

Main results for the dispersion and the spectral func- 
tions are not changed much in comparison with the pre- 
vious ones as shown in Fig. 1131 Larger hybridization 
between the subbands at small value of U e ff results in 
increase of the dispersion and the intensity of the up- 
per one-hole subband. This trend is also seen in the 
DOS in Fig. HH At weak doping the Mott gap be- 
tween the subbands is observed despite the intermediate 
Coulomb energy U e ff = At, only a half of the bare band- 
width W ~ 8*. This can be explained by a reduction 
of the bandwidth caused by strong spin correlations in 
the undcrdoped region up to W ~ 8\t'\ as discussed in 
Sect. Ill Bl below the equation (|15fl . In the overdoped 
case at d = 0.3 when the spin correlations become weak 
the gap between the subbands vanishes. 

Noticeable changes are observed for the FS shown in 
Fig. 1151 and in Fig.^3 In the first plot where the FS was 
determined by the equation (f^Uf) we see a large pocket 
at small doping 5 = 0.1 which opens with doping or 
temperature increase. At the ovcrdoping for 5 = 0.3, 
the FS transforms to the electron-like as in the previous 
calculations. This transformation is confirmed by cal- 
culations of the electron occupation numbers shown in 
Fig. El It should be noted that a pronounced hole 
pocket in the new set of the model parameters is caused 
by the t" contribution which results in a large dispersion 
in the (n, 0) — ► (0, n) direction (oc t" (cos2fca; + cos2fc y )) 
disregarded in the previous set of the parameters. A re- 
markable feature of these results is that the part of the 
FS close to the T(0, 0) point in the nodal direction in 
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FIG. 13: (Color online) Dispersion curves for A = 4t along the 
symmetry directions T(0, 0) -> M(w,tt) -> X(tt,0) -> T(0,0) 
at S = 0.05 (upper panel) and 8 — 0.3 (bottom panel). 
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FIG. 14: (Color online) Doping dependence of the DOS for 
A = it. 



Fig. 1151 does not shift much with doping (or tempera- 
ture) being pinned to a large FS as observed in ARPES 
experiments (see, e.giSS). In fact, only this part of the FS 
was detected in ARPES experiments where the spectral 
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FIG. 15: (Color online) Doping dependence of the FS for 8 = 
0.1 (full line at T = 0.0M and dotted line at T = 0.3i), 8 = 0.2 
(dashed line), and <5 = 0.3 (dot-dashed line) for A = At. 
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FIG. 16: (Color online) A(k,w = 0) on the FS at 8 = 0.05 
(left panel) and 8 = 0.1 (right panel) at T = 0.03i for A = At. 

function A e z(k, u> = 0) shown in Fig. Ilfil was measured. 

Concerning the self-energy effects and kinks, they are 
similar to the case for A = 8i and confirm a strong 
influence of spin correlations on the QP spectra renor- 
malization. As shown in Fig. the coupling constant 
A = — (<9Re£(k, uo)/duj) u= o being large at small dop- 
ing distinctly decreases with ovcrdoping at 5 — 0.3 ac- 
companied by suppression of the imaginary part of the 
self-energy. In conclusion, the alternative set of parame- 
ters with a moderate effective Coulomb energy U e f f = At 
in the Hubbard model Q confirms an important role of 
AF correlations in the electronic structure of system with 
large single-site Coulomb interaction. 



IV. CONCLUSION 

In the present paper the theory of electronic spectra 
in the strong correlation limit for the Hubbard model 
in a paramagnetic state has been formulated. By em- 
ploying the Mori-type projection technique for the ther- 
modynamic GFs in terms of the Hubbard operators, 
we consistently took into account charge carrier scatter- 
ing by dynamical spin fluctuations and derived the self- 




FIG. 17: (Color online) The electronic occupation numbers 
iV k at T = 0.03t for 8 = 0.05 (upper panel) and at 8 = 0.3 
(bottom panel) for A = At. 



consistent system of equations for the GF (|23() and the 
self-energy <|22|) evaluated in the NCA which neglects the 
vertex corrections. Though in the Hubbard model l|T)l 
the electron coupling to spin-fluctuations is not weak, 
it is of the order of the hopping parameter, the vertex 
corrections should not be so important in this case due 
to kinematic restrictions imposed on the spin-fluctuation 
scattering. As was shown for the t-J model^, the lead- 
ing two-loop crossing diagram identically vanishes, while 
the next three-loop crossing diagram gives a small contri- 
bution to the self-energy. In any case, the NCA for the 
self-energy can be considered as a starting approxima- 
tion for a model with strong coupling. As we discussed 
at the end of Sect. IIII Al the self-consistent system of 
equations for the self-energy in the classical limit in our 
approach are similar to the two-particle self-consistent 
approach (TPSC)2£ or the model of short-range static 
spin (charge) fluctuationsa^. Numerical results for the 
spectral density and the FS in the NCA approximation 
for the self-energy are quite similar to the studies within 
the generalized DMFT^i^ where all diagrams for elec- 
tron scattering by spin (or charge) fluctuations in the 
static approximation were taken into account. Our re- 
sults are also in accord with calculations based on the 
cluster approximation^ and the TPSG2£. 

In the present paper we have not presented a fully self- 
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FIG. 18: (Color online) Energy dependence of the real and 
imaginary parts of the self-energy E(k, uj) for A = 4t at the 
r(0,0), S(tt/2, 7r/2) and M(tt,it) points at S = 0.1 (upper 
panel) and S = 0.3 (bottom panel). 



consistent theory for the single-electron GF and the dy- 
namical spin and charge susceptibility. This demands 
rather involved calculations of the collective spin and 
charge excitation spectra which is beyond the scope of 
the present paper. Instead, we have used a model for the 
dynamical spin susceptibility (|26[) which is usually cm- 
ployed in phenomcnological approach. However, a varia- 
tion of the electron (hole) interaction with spin fluctua- 
tions in our theory is strongly restricted since the vertex 
of the interaction is given by the hopping parameters @ 
in the Hubbard model, while an intensity of spin fluctu- 
ations at the AF wave- vector Q (C(£) in the Table UJ) 
determined by the AF correlation length £ is fixed by the 
sum rule l|27|) . A variation of the cut-off energy cu s does 
not affect noticeably the numerical results, as we have 
checked. The resulting coupling constant A obtained in 
our calculations (see Sect. IIII D|) seems to be too large 
in comparison with ARPES results. This discrepancy 
can be caused by disregarding scattering on charge fluc- 
tuations in the dynamical susceptibility model (|23|) and 
electron-phonon interaction which may reduce the con- 
tribution from the electron-spin interaction. 

The main conclusion of the present study is that a de- 
cisive role in renormalization of the electronic spectrum 
in strongly correlated system as cuprate superconductors 
is played by electron interaction with spin-fluctuations 
which is in accord with other studies (e.g.^iSii^). The 
numerical results for the electron dispersion in Sect. IIII Bl 
the FS and the occupation numbers in Sect. IIII HI and the 
self-energy in Sect. IIII Dl unambiguously approved this 
conclusion. With doping or temperature increasing, spin 
correlations are suppressed which results in transition 
from a strong to a weak correlation limit. These observa- 
tions were confirmed also by a consideration of the model 
with intermediate Coulomb correlations in Sect. IIII El 

A theory of superconducting transition within the 
present theory will be considered elsewhere. 
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